Using the cluster mass function from weak lensing to constrain neutrino masses 
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We discuss the variation of cosmological upper bounds on Mi,, the sum of the neutrino masses, 
with the choice of data sets included in the analysis, pointing out a few oddities not easily seen when 
all data sets are combined. For example, the effect of applying different priors varies significantly 
depending on whether we use the power spectrum from the 2dFGRS or SDSS galaxy survey. A 
conservative neutrino mass limit of Mi, < 1.43eV (95%C.L.) is obtained by combining the WMAP 
3 year data with the cluster mass function measured by weak gravitational lensing. This limit has 
the virtue of not making any assumptions about the bias of luminous matter with respect to the 
dark matter, and is in this sense (and this sense only) bias-free. 



I. INTRODUCTION 

The fact that neutrinos undergo flavour oscillations, 
implying that not all neutrinos are massless, is one of 
the most important discoveries in particle physics in the 
last decade. Oscillations are only sensitive to the mass- 
squared differences between the neutrino mass eigen- 
states (see e.g. [1113])' although these are suggestive 
of the overall mass scale, they cannot pin the absolute 
mass scale firmly. For that, one needs probes like tri- 
tium beta decay [3] or neutrinoless double beta decay . 
At the moment, however, the strongest upper bound on 
the neutrino mass scale comes from cosmology (see e.g. 
[i, i]). The cosmolo gical mass limits are based on the 
implications of neutrino masses for clustering of matter, 
in particular the fact that neutrinos can free-stream out 
of density perturbations and impede structure formation 
on small scales if they are a significant fraction of the 
dark matter. For a thorough review, see Q. 

Impressive as they are, the cosmological neutrino mass 
limits involve several assumptions. First of all, they as- 
sume that the underlying cosmological model is of the 
Friedmann-Robertson- Walker type, that gravity is de- 
scribed by general relativity, the primordial power spec- 
trum of density perturbations is a scale-free power-law, 
and that the dark energy is a cosmological constant. 
Modifications of this simple picture have been consid- 
ered 0, and in particular the degeneracy between neu- 
trino masses and the dark energy equation of state has 
been investigated in some recent papers d, [l^, El, ■ 
Secondly, the upper limit depends on the cosmological 
data sets used in the analysis. It is this last aspect of the 
problem we will consider in the present paper, working 
within the context of the standard ACDM paradigm. 

The angular power spectrum of the temperature 
anisotropics in the cosmic microwave background (CMB) 
radiation is generally considered to be the cleanest cos- 
mological data set. Unfortunately the CMB temperature 
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power spectrum cannot provide upper limits on bet- 
ter than ~ 1.6 eV [l3|, which is almost reached already 
now with the WMAP 3-year data [l3,[ll[ll. 

An im- 
provement of this limit requires some probe of the matter 
distribution. The most common probe is the power spec- 
trum of the galaxy distribution, as determined by large 
redshift surveys like the 2 degree Field Galaxy Redshift 
Survey (2dFGRS) [iT^ and the Sloan Digital Sky Survey 
(SDSS) [13] ■ The distribution of the galaxies is biased 
with respect to that of the dark matter, but if the bias 
is independent of length scale, this is not a serious limi- 
tation. However, recent results suggest that the story of 
bias might be more complicated than the simple assump- 
tion of scale-independent bias that has gone into earlier 
cosmological mass limits [l3, [13] • 

Ideally, one would like to rely on probes that are sen- 
sitive to the total matter distribution. One such probe 
is measurements of weak gravitational lensing by mat- 
ter along random lines of sight [lli[23]. However, even 
the most ambitious of these 'cosmic shear' surveys have 
covered less than one percent of the full sky, thus prob- 
ing only a limited cosmological volume. A somewhat 
different approach is to count the number of rare, high 
density peaks in the matter fiuctuations within a larger 
volume. These peaks correspond to very massive clus- 
ters of galaxies, their abundances (the so-called cluster 
mass function) being highly sensitive to the amplitude of 
the matter power spectrum (see e.g. [23| and references 
therein) on corresponding scales. Until recently, cluster 
masses were generally derived based on observables of the 
baryonic mass component, such as the temperature of the 
X-ray emitting intra-cluster medium (ICM) . In this case, 
one had to calibrate the X-ray temperature-mass rela- 
tionship, either from simulations or from observational 
mass measurements relying on assumptions about hy- 
drostatical equilibrium of the ICM. However, Dahle [2J] 
provided for the first time a measurement of the clus- 
ter mass function (CMF) where the cluster masses have 
been established directly by weak gravitational lensing. 
The main result of this paper is an upper bound on the 
sum of the neutrino masses obtained by combining the 
WMAP 3-year data with the CMF derived by Dahle [H]. 
Although not as impressive as e.g. the limit obtained 
by combining all available cosmological observations, in- 
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eluding the Lyman alpha forest power spectrum (see e.g. 
BUEl), this limit is robust in the sense that it makes use 
of clean cosmological probes, and involves a minimum of 
assumptions. 

II. MASSIVE NEUTRINOS IN COSMOLOGY 

We work within the standard cosmological paradigm 
of the ACDM model, and leave studies of significant de- 
viations from this model for a future study. Our adopted 
notation is as follows: We denote as ili the present den- 
sity of component i in units of the density of a spatially 
flat universe. Note that flm is the total density of all non- 
relativistic components, so that Q^a = ^^cdm + + ^i^, 
with fib being the present baryon density. The density 
parameter fl^ is determined by the sum of the neutrino 
masses and the present value of the Hubble parame- 
ter Ho = lOO/i kms-i Mpc"^ [1^: 



As we will see, present cosmological probes (with the 
possible exception of the Lyman alpha forest Q ) are only 
sensitive to neutrino masses larger than a few tenths of 
an eV. In this regime, the neutrino mass spectrum would 
be degenerate, and hence = 3m^, where rrii, is the 
common mass of the neutrino mass eigenstates. 

Within the ACDM model the effect of massive neu- 
trinos on structure formation is simple. Their influ- 
ence on structure formation is governed mainly by the 
quantity — il^/fim. While relativistic, the neutri- 
nos free-stream out of density perturbations on scales 
smaller than the Hubble radius. Once they become non- 
relativistic, the free-streaming scale is set by their root- 
mean-square velocity. The maximum length scale below 
which all scales are affected by neutrino free-streaming is 
thus set by the horizon size when the neutrinos became 
non-relativistic. This quantity is given by Ji] 

^ 0.018 [^y^'nli'hMpc-' 

= OAfy^fln^hkMpc-^. (2) 

Note that the combination f2,„/i enters this expression; 
this quantity determines the Hubble radius at matter- 
radiation equality and is an important length scale in 
the power spectrum of matter density fluctuations. This 
scale is similar to the neutrino free-streaming scale: be- 
low this scale, density perturbations are suppressed. 
Hence, there is a degeneracy between fl^h and 

For scales below the neutrino free-streaming scale, 
k > the main effect is a reduction by a factor 

1 — of the source term in the equation for the lin- 
ear growth of density perturbation. Hence, the amount 
of suppression is also set by the parameter /y, and so it 
is this, and not Mi, which is directly accessible in mea- 
surements of large-scale structure. Since we can write 




klhMpc'l 



FIG. 1: Illustration of how combining data sets break de- 
generacies. The left panel shows the matter power spectrum 
in four different models, where the baryon density has been 
artificially set low and a scale invariant power spectrum of adi- 
abatic primordial perturbations has been assumed. The black 
lines are for models with Q^, = 0. The full line has Sim = 0.25 
and h = 0.735, whereas the dotted line has flm ~ 0.50 and 
h — 0.3675. The red lines are models with massive neutrinos. 
The dashed line has fim = 0.25, fii/ = 0.05 and h — 0.735, and 
the dot-dashed fine has Sim ~ 0.50, fi^ — 0.10 and h = 0.3675. 
Note that = 0.2 in both cases. The right panel shows the 
CMB power spectra for the two models with massive neutri- 
nos, with the same line coding. The models have very dif- 
ferent CMB power spectra, and hence measurements of the 
CMB anisotropies help constraining Q^. 

Mj, = 93.14/^r2ni/i^ eV, we see that one also needs a 
constraint on ^j^h^ to go from a limit on f^, to a limit on 
Mj^. Figure [1] illustrates this situation in a toy example. 
It is seen that the matter power spectrum must be com- 
bined with other probes, for example the CMB power 
spectrum, in order to provide stringent constraints on 



III. DATA AND METHODS 

In our analysis we have used data from both CMB, 
observations of large scale structures (LSS), type la su- 
pernovae (SNIa), baryon acoustic oscillations (BAO), 
and additional priors on the Hubble parameter and the 
baryon content of the universe. We have also applied 
constraints on the cluster mass function from weak grav- 
itational lensing. 

A. CMB data 

The 3-year data release from the WMAP team is at 
present the most constraining set of CMB observations. 
The WMAP experiment is a full-sky survey of CMB tem- 
perature anisotropies. With this release, a Fortran 90 
code for calculating the likelihood of a given CMB power 
spectrum against the data was also provided. However, 
in ref. [2^ it was noted that this likelihood code applied 
a sub-optimal likelihood approximation on large scales 
{I ^ 30), and in ref. [1^ the authors pointed out an over- 
subtraction of unresolved point sources on small angular 
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scales (/ > 30). Correcting for these two effects resulted 
in a slight shift in some cosmological parameter values. 
In ref. it was shown that these corrections to the 

WMAP likelihood code also shifted the upper limit on 
the neutrino masses ft'om Mi, < 1.90eV to < 1.57eV 
when using WMAP data only. In a revised version of 
the WMAP likelihood code 47j, these corrections have to 
some extent been accounted for, and the results using 
this code is in better agreement with refs. [l^, than 
the previous version. In this paper we will make use of 
this revised likelihood code from the WMAP team. We 
note, however, that the upper limits on Mi, would have 
been slightly lower with the code used in refs. (TgI. [29j. 

In this work we are not considering CMB data from 
other experiments than WMAP, since we want to restrict 
the number of data sets as much as possible, and since 
these additional data sets have been shown not to im- 
prove the limits on Mi, |l6| . 



B. Large scale structure 

Large scale structure surveys probe the matter distri- 
bution in the universe by measuring the galaxy-galaxy 
power spectrum Pg{k,z) = (|5g(fc, Since mas- 

sive neutrinos give a distinct imprint on the LSS due 
to their free-streaming effect, data from galaxy surveys 
have proven to be important to put tight constraints on 
neutrino masses. However, such surveys are also troubled 
by bias effects that are hard to quantify. For example, 
it is assumed that in the linear perturbation regime, the 
total matter power spectrum, Pm, is proportional to the 
galaxy power spectrum by the simple relation Pg = b^Pm- 
Recently, the use of this simple relation has been debated 
[20j , and it is at present unclear what the precise cor- 
rections from a scale-dependent b might be. 

There are two galaxy surveys of comparable size, the 
2 degree Field Galaxy Redshift Survey (2dFGRS) [13, 
and the Sloan Digital Sky Survey (SDSS) main sample 
[30t . As has been pointed out in e.g. [H, [U, H^l , there 
seems to be some tension between the results from 2dF- 
GRS and SDSS. As a relevant example for this paper, the 
combination SDSS-I-WMAP prefers larger values of both 
f2in and trg (the rms mass fluctuations in spheres of ra- 
dius 8ft.~^ Mpc) than what is found when using data from 
2dFGRS-f WMAP. In a recent analysis of the 2dFGRS- 
SDSS discrepancy [34I, the authors conclude that it is 
caused by non-linear bias effects. They also claim that 



SDSS is more sensitive to these effects than the 2dFGRS 
sample, since the SDSS sample contains more red galax- 
ies, which are believed to cluster slightly different from 
blue galaxies. 

Regardless of their origin, the reported inconsistencies 
between the 2dFGRS and SDSS galaxy surveys tell us 
that we should always be cautious when including LSS 
data in our parameter estimation, at least until the scale 
dependence of the bias parameter is better understood. 

When including 2dFGRS and SDSS data in our analy- 



sis we have tried to use as commonly used analysis tech- 
niques as possible. We have therefore stuck to the default 
analysis and parameter choices in the CosmoMC code. 
For the 2dFGRS data this implies discarding all scales 
with k > 0.15/iMpc~^. This corresponds to the choice 
made in e.g. [3l|. For the 2dFGRS data we also use 
the prescribed correction for nonlinearities given in p^ . 
Pg = (1-f Qfc2)/(l-|-^fc) with A = 1.4 and Q = 4.0 ±1.5. 
For SDSS we discard scales with k > 0.20/iMpc~^, as 
is done in e.g. [l^. In this paper they find that using 
this cutoff scale yields results that are consistent with 
the what they find when applying more a more conser- 
vative cutoff on k. No corrections for non-linearities are 
included in our SDSS analysis. The treatment of cut- 
offs and nonlinearities in the analysis of the 2dFGRS 
and SDSS data are different. Thus, if any tension is 
found between these two data sets here, it is not obvious 
whether this stems from problems within the datasets or 
if it is caused by the way the data sets are analysed here. 
However, any such tension would illustrate that there are 
problems related to using LSS surveys for parameter es- 
timation when adapting these commonly used analysis 
methods for the two surveys. 



C. Type la supernovae 

The luminosity distance-redshift relationship mea- 
sured by observations of supernovae of type la (SNIa) 
provides the most direct evidence for cosmic acceleration 
and dark energy. There are still open questions regarding 
both the exact mechanism behind these supernovae and 
their use as 'standard candles', and so we choose not to 
rely on SNIa in our robust neutrino mass limit. When we 
do use them, we use the data from the Supernova Legacy 
Survey (SNLS) [33,]. 



D. The cluster mass function 

Massive galaxy clusters are extremely rare high-density 
peaks in the matter fluctuations, containing matter orig- 
inating from a co-moving volume spanning ^ 10h~^ Mpc 
which has undergone gravitational collapse. Their abun- 
dances, as quantified by the CMF, are thus a sensitive 
probe of erg and fJm- While all previous measurements 
of the CMF have been based on observations of baryonic 
tracers of cluster mass (with inherent uncertainties and 
possible biases which are not yet fully understood; see 
e.g. [131), the CMF of [l^l was derived from weak gravi- 
tational lensing measurements of the masses of a volume- 
limited sample of X-ray luminous clusters. The details of 
the data reduction and weak gravitational shear estima- 
tor are given in and [13] describes the derivation of 
cluster masses from gravitational lensing measurements. 
The cluster sample of was selected from a large vol- 
ume of 8.0 X 10^(/i~^Mpc)'^ (assuming a spatially flat 
universe with = 0.3), above a threshold value in X- 
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ray luminosity (ix)- Using Lx as a proxy for mass in the 
cluster selection could in principle introduce a baryonic 
bias to the CMF measurements, depending on Lx mea- 
surement uncertainties and the intrinsic scatter around 
the mass-Lx relation. In [24J, any such effects were ef- 
fectively removed, by carefully calibrating the amplitude 
of the mass-Lx relation and the scatter around the mean 
relation (thereby estimating the sample completeness as 
a function of mass), and by calculating cosmological pa- 
rameter constraints exclusively based on the clusters well 
above the mass threshold corresponding to the X-ray lu- 
minosity threshold of the sample (where the sample is 
virtually complete) . Other statistical and systematic un- 
certainties in the derived cluster masses and the effects of 
these uncertainties on the CMF constraints are discussed 
in detail in [23 |. 

Parameter constraints in the (Tg — flm plane were ob- 
tained by fitting the observed cluster abundances in three 
mass intervals to theoretical predictions for the CMF 36] . 

We have found that the value of from the CMF from 
[2^ can be approximated by the fit-function 

XcMF = 10000u'^-I-6726m3+1230u2-4.09u-I-0.004, (3) 

where u = CT8(ilin/0.3)'' '^^ — 0.67. This result is rather 
insensitive to the choice of theoretical mass function: As 
noted by |2j|, replacing the CMF prediction used here 
[Sfii with a more recent predicted CMF based on N-body 
simulations [33| (taking into account the differing mass 
definitions in these two works) only results in a shift of 
0.01 in u. 



where Dm{z) is the comoving angular diameter distance, 
the authors impose the constraint 

(jl \ —0.35 
^) (1 + 0.94/,) ±0.017. (5) 

Throughout our analysis we have also applied a top-hat 
prior on the age of the universe, lOGyr < Age < 20Gyr. 

F. Parameter estimation 

For the parameter estimations we have used 
the publicly available Markov chain Monte Carlo 
code CosmoMC [i^. We have used a basic 
seven-parameter model, varying the parameters 
{nbh'^,n^,\og[W^As),h,ns,T,M^). is defined in 
eq. ([1]). For the other parameters, the exact definitions 
are given by the CosmoMC code. We have assumed 
spatial flatness, no running of the scalar spectral index, 
that the dark energy is a cosmological constant and 
that the tensor to scalar fluctuation amplitude is negli- 
gible. These assumptions are well motivated by current 
available data [l3|- Also, adding these extra degrees of 
freedom do not affect the limits on M^, drastically 1^. 

Often cosmological neutrino mass limits are found us- 
ing a large number of different data sets and priors si- 
multaneously. In this analysis we consider in total 48 
different combinations of data sets and priors to see how 
the different data sets alter the neutrino mass limit. 



IV. RESULTS 



E. Other priors 

We have also tested the sensitivity of neutrino mass 
limits on priors from the Hubble Space Telescope Key 
Project on the Extragalactic Distance Scale (HST), and 
Big Bang nucleosyntesis (BBN) predictions. 

From HST, we use a prior on the Hubble parameter of 
h = 0.72 ± 0.08 [Hi, and from BBN we use a prior on 
the physical baryon density today, Q.bh? — 0.022 ± 0.002 
(39. . ,40, . ,41i] . There are hints of some tension between 
the WMAP constraint on Q.bh? and the value of this 
quantity inferred from the "'He abundance [13], so in our 
robust limit we will not use the BBN prior. Since there is 
still some debate about the value of the Hubble constant 
derived from the HST Key Project [H, SI], we also drop 
this prior when deriving our robust limit. 

We have also added information on the position of the 
baryonic acoustic oscillation (BAO) peak in the luminous 
red galaxy (LRG) sample in the SDSS survey [i^. This 
prior is implemented by an effective flt function as de- 
scribed in [5|]. Using an effective parameter 



A 



BAO 



DM{.zf 



H{z) 



1/3 



Hi 



(4) 



The neutrino mass limits found in our analysis are 
summarized in Table [D The limits quoted in the ta- 
ble range from < 6eV from using LSS data only, 
via < 1.75eV from WMAP data only, and down to 
AU < 0.40eV for a combination of WMAP, SDSS, SNLS 
BAO and HST data. 



A. WMAP + priors 

As one might have expected, the inclusion of WMAP 
data is crucial for finding a good neutrino mass limit. The 
upper limit from WMAP data alone found here. Mi, < 
1.75eV, resides between the former results found by the 
WMAP team for their 3 year data, < 2.0eV [l^. 
and the results from ref. 16"| using the same data with a 
modified likelihood code, < 1.57eV. In this case we 
see that both a CMF and HST prior will help significantly 
in constraining M^, and adding the CMF prior to the 
WMAP data yields an upper bound of AU < 1.43eV. 

The CMF prior put constraints on a combination of 
rim and (Ts J and in Figure [2] we show confidence contours 
in the plane of M^, and cr8(f^m/0.3)°-37 from WMAP data 
only, and when we add the CMF prior. It is well known 
that there is a positive correlation between Al^, and 
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Data set 



no prior CMF HST BBN 



WMAP 1.75eV 

WMAP+2dFGRS 1.02eV 

WMAP+SDSS 1.05eV 

WMAP+SNLS l.lOeV 

SDSS 5.8eV 

2dFGRS 5.2eV 

WMAP+2dFGRS+SNLS O.SleV 

WMAP+SDSS+SNLS 0.44eV 

WMAP+2dFGRS+SNLS+BAO 0.84eV 

WMAP+SDSS+SNLS+BAO 0.42eV 

WMAP+SDSS+BAO 0.55eV 

WMAP+2dFGRS+BAO 1.04eV 



1.43eV 1.47eV 1.73eV 
0.89eV 0.92eV l.OleV 
1.13eV 0.87eV 1.05eV 
0.70eV 1.02eV 1.07eV 
4.6eV e.leV 6.0eV 
5.3eV 5.3eV 5.2eV 
0.64eV 0.76eV 0.80eV 
0.57eV 0.42eV 0.44eV 
0.69eV 0.79eV 0.84eV 
0.56eV 0.40eV 0.42eV 
0.76eV 0.52eV 0.56eV 
0.89eV 0.92eV l.OleVFIG. 
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2: 68% and 95% C.L. contours in the plane of 



TABLE L Estimated 95% C.L. upper limits on M^. The 
CMF, HST and BBN priors are added one at the time to the 
combination of data sets given in the left column. 



My and the effective parameter from the CMF constraint, 
a8(f2m/0.3)'''^'^. The dashed black line shows the contours 
when using WMAP data only, and the solid black lines show 
the results when adding the prior from the CMF. 



in the CMB power spectrum. Both these parameters 
will alter the amplitudes of the CMB peaks, by shift- 
ing the time of matter-radiation equality. At this time 
neutrinos in this mass range were still relativistic, and 
contributed to the radiation part. Keeping constant 
and increasing will thus postpone the time of matter- 
radiation equality and therefore enhance the amplitude 
of the acoustic peaks. To shift the time of equality back 
and lower the acoustic peaks, one has to increase rim cor- 
respondingly. However, the effective parameter from the 
CMF constraint, CT8(r2m/0.3)'^ '^^, also contains the am- 
plitude parameter erg which is negatively correlated with 
Mu (from the same reasoning) . From Figure [5] we see 
that this negative correlation is strong enough to make 
also the correlation between and crgC^^m/O.S)" '^^ neg- 
ative; that is, small values of (T8(ilin/0.3)'' '^^ favor a large 
neutrino mass. When we add the CMF prior the allowed 
region of low values of (T8(f^m/0.3)°-^'' shrinks, and the 
upper limit on Mn is reduced correspondingly. 

Note that with this limit (Aljj < 1.43eV) we have not 
used any LSS or SNIa data, or any other prior on Hq or 
ilb- We therefore claim this neutrino mass limit to be 
robust. 



B. Including large scale structure 

When adding LSS data the neutrino mass limit 
improves significantly, giving Af^ < 1.02eV for 
WMAP-|-2dFGRS and < l.OSeV for WMAP -f SDSS. 
However, when adding the CMF prior, a strange effect 
occurs. While this improves the M^, limit to < 
0.89eV in the case of WMAP+2dFGRS, the Hmit in- 
creases to < l.lSeV with WMAP-hSDSS. This may 
indicate a inconsistency between the CMF prior and 
our SDSS analysis. In Figure [3] we show the marginal- 
ized distribution of the fT8(f^m/0.3)°'^^ parameter 
when using WMAP, WMAP-I-CMF, WMAP+2dFGRS, 



WMAP+SNLS or WMAP-^SDSS. For the latter, the dis- 
tribution deviates significantly from the four former, and 
we see that the WMAP-I-SDSS data set prefers larger 
values of as{^m/0.3)'^'^'' than than the other combina- 
tions. This is in good accordance with the results ob- 
tained in ref. [3l|, where they found that WMAP-hSDSS 
preferred larger values for both fl^ and ag than what was 
found with WMAP alone or WMAP+2dFGRS (using the 
WMAP 1-year data) . Figure [3] indicates that we should 
be careful when using data from SDSS in combination 
with 2dFGRS or CMF. Doing this, the inconsistencies 
may produce artificially narrow parameter distributions. 
In our case, the result of using WMAP+SDSS+CMF is 
that we not only get an upper limit on AI^, but also 
a lower limit, such that our 95% C.L. limit becomes 
0.23eV < Ml, < 1.13eV for this combination of data 
sets. In Figure [3] we show the confidence contours in the 
plane of asi^njO-^f -^'' and for WMAP+2dFGRS 
and WMAP-I-SDSS with and without the CMF prior. 
From these plots it is clear how the incompatibility of 
SDSS and CMF produces the lower limit on M^. 

When adding LSS data to the WMAP data set, we see 
that the HST prior is still important to improve the M^, 
limits, while the BBN prior is not needed. 

Also, it is interesting to note that adding the BAO 
prior to the WMAP-|-2dFGRS data set combination has 
essentially no effect for the neutrino mass limits, whereas 
adding the BAO prior to the WMAP -f SDSS data sets 
improves the neutrino mass limit by almost a factor two. 
Again this illustrates the differences between the 2dF- 
GRS and SDSS sample in our analysis. The reason why 
the BAO prior is more important when using SDSS, is 
the larger preferred value of fim/i^ for the SDSS sample. 
Adding the BAO prior will constrain the allowed region 
of large ilm^^, and thus also the allowed space for large 

We have also tried to constrain M^, using LSS data 



6 



WMAP 




0.4 0.5 0.6 0.7 „„ 0.8 0.9 1 



FIG. 3: Marginalized one-dimensional distribution of 
(78 (r2m/0.3)°'^'' when using WMAP data (thick black line), 
WMAP+CMF (solid red line), WMAP+2dFGRS (dotted 
blue line), WMAP+SDSS (dash-dotted green line) and 
WMAP-FSNLS (dashed purple line). 




FIG. 4: 68% and 95% C.L. contours in the plane 
of asiflm/0-3)°-^'' and M^. Left panel: Con- 

tours from WMAP-I-SDSS (black, dashed line), and 
WMAP-HSDSS-HCMF (red, solid line). Right panel: Con- 
tours from WMAP-|-2dFGRS (black, dashed line), and 
WMAP4-2dFGRS-hCMF (red, solid line). 

alone, resulting in mass limits of order AI^, ^ 6eV. So 
although LSS data in principle is a sensitive probe for 
neutrino masses, one needs to add CMB data to break 
parameter degeneracies. 



C. Including supernova data 

The inclusion of supernova data turns out to be as im- 
portant as LSS for constraining neutrino masses. Again, 
this can be understood by the M,y-r2i„ degeneracy. SNIa 
data is an effective probe of the amount of dark energy in 
the universe, and under the flatness assumption this will 
also automatically constrain U^. In Figure [3] we see that 
also the data set combination WMAP+SNLS seems to 
be consistent with WMAP alone, WMAP-|-2dFGRS and 
CMF prior, whereas it seems to be in some tension with 



our SDSS analysis. In the case of WMAP-f SNLS we see 
from Table U that adding the CMF prior improves the 
neutrino mass limit significantly, and that the inclusion 
of the CMF prior is a lot more constraining for neutrino 
masses than both the HST and BBN prior in this case. 
This means that if one believes in the SNIa measure- 
ments, we have an upper limit of M^, < 0.70eV without 
using LSS measurements at all. 

Note also that we can get a neutrino mass limit as low 
as Ml, < 0.40eV without using Ly-a data, if we com- 
bine WMAP+SNLS+SDSS+BAO+HST. But this limit 
weakens to < 0.79eV by substituting SDSS with 2dF- 
GRS. Again, this illustrates how sensitive the neutrino 
mass limit is to the choice of LSS sample and combina- 
tion of priors. 



V. CONCLUSION 

In this paper we have studied cosmological neutrino 
mass limits, and how sensitive these limits are to different 
choices of data sets and priors. We have also included a 
new prior from the cluster mass function measured by 
weak gravitational lensing, which is a direct probe of the 
total mass distribution in the universe. 

We report neutrino mass limits using 48 different com- 
binations of data sets and priors, and we find that the 
neutrino mass limit is very sensitive to small changes in 
combinations of data sets and priors. Especially striking 
results are found when interchanging data sets between 
the SDSS and 2dFGRS galaxy surveys. For example 
will the combination WMAP-hSDSS-hBAO give A-U < 
0.52eV at 95% C.L., while WMAP-t-2dFGRS+BA0 gives 
My < 0.97eV. These discrepancies occur because of a 
slight inconsistency between the SDSS data as analysed 
here and many of the other data sets used. Combining 
inconsistent data sets may of course lead to unreliable re- 
sults. Also, the tension between the two galaxy surveys 
indicate that there may be systematics related to e.g. the 
scale dependence of the bias parameter that have to be 
better understood before we can rely fully on the LSS 
analysis. By discarding these data sets, and in addition 
refraining from using constraints from SNIa, HST, and 
BBN, we end up with a conservative, although robust 
cosmological neutrino mass limit of M^, < 1.43eV from 
using only WMAP data in combination with the CMF 
prior. 



Acknowledgments 

0E, JRK and HD acknowledge support from the 
Research Council of Norway through project numbers 
159637, 162830 and 165491. Some of the results in this 
work are based on observations made with the Nordic 
Optical Telescope, operated on the island of La Palma 
jointly by Denmark, Finland, Iceland, Norway, and Swe- 



7 



den, in the Spanish Observatorio del Roque de los Mucha- chos of the Institute de Astrofisica de Canarias. 



[1] M. Maltoni, T. Schwctz, M. A. Tortola, and J. W. F. 

Valle, New J. Phys. 6, 122 (2004), hcp-ph/0405172. [26] 
G. L. Fogli et al. (2006), hcp-ph/0608060. 

K. Eitel, Nucl. Phys. Proc. Suppl. 143, 197 (2005). [27 
S. R. Elliott and P. Vogel, Ann. Rev. Nucl. Part. Sci. 52, [28^ 
115 (2002), hep-ph/0202264. [29 
A. Goobar, S. Haimestad, E. MortscU, and H. Tu (2006), 
astro-ph/0602155. [30 
U. Seljak, A. Slosar, and P. McDonald (2006), astro- 
ph/0604335. [31 
J. Lesgourgues and S. Pastor, Phys. Rept. 429, 307 
(2006), astro- ph/0603494. [32 

C. Zunckcl and P. G. Forrcira (2006), astro-ph/0610597. 
S. Hannestad, Phys. Rev. Lett. 95, 221301 (2005), astro- [33 
ph/0505551. 

K. Ichikawa and T. Takahashi (2005), astro-ph/0510849. [34 
A. De La Macorra, A. Melchiorri, P. Serra, and R. Bean 
(2006), astro-ph/0608351. 

J.-Q. Xia, G.-B. Zhao, and X. Zhang (2006), astro- 
ph/0609463. [35 
K. Ichikawa, M. Fukugita, and M. Kawasaki, Phys. Rev. 
D71, 043001 (2005), astro-ph/0409768. [36 

D. N. Spcrgcl et al. (2006), astro-ph/0603449. 
M. Fukugita, K. Ichikawa, M. Kawasaki, and O. Lahav, [37 
Phys. Rev. D74, 027302 (2006), astro-ph/0605362. 
J. R. Kristiansen, H. K. Eriksen, and O. Elgaroy (2006), [38 
astro-ph/0608017. 

S. Cole et al. (The 2dFGRS), Mon. Not. Roy. Astron. [39 
Soc. 362, 505 (2005), astro-ph/0501174. 

M. Tegmark et al. (SDSS), Phys. Rev. D69, 103501 [40 

(2004), astro-ph/0310723. 

W. J. Percival et al. (2006), astro-ph/0608636. [41 
R. E. Smith, R. Scoccimarro, and R. K. Shcth (2006), 
astro-ph/0609547. [42 
A. Refregier, Annual Review of Astronomy & Astro- 
physics 41, 645 (2003), astro-ph/0307212. [43 
P. Schneider, in Saas-Fee Advanced Course 33: Grav- 
itational Lensing: Strong, Weak and Micro, edited by 
G. Meylan, P. Jetzer, P. North, P. Schneider, C. S. [44 
Kochanek, and J. Wambsganss (2006), pp. 269-451. [45 

G. M. Voit, Reviews of Modern Physics 77, 207 (2005), 
astro-ph/0410173. [46 

H. Dahle Astrophys. J. 653, 954 (2006), astro- 
ph/0608480. [47 
M. Cirelh and A. Strumia, JCAP 0612, 013 (2006), 



astro-ph/0607086. 

G. Mangano et al., Nucl. Phys. B729, 221 (2005), hep- 
ph/0506164. 

G. Hinshawet al. (2006), astro-ph/0603451. 

H. K. Eriksen et al. (2006), astro-ph/0606088. 

K. M. Huffenbergcr, H. K. Eriksen, and F. K. Hansen 

(2006), astro-ph/0606538. 

M. Tegmark ct al. (SDSS), Astrophys. J. 606, 702 (2004), 

astro-ph/0310725. 

A. G. Sanchez ct al., Mon. Not. Roy. Astron. Soc. 366, 
189 (2006), astro-ph/0507583. 

S. Cole, A. G. Sanchez, and S. Wilkins (2006), astro- 
ph/0611178. 

P. Asticr ct al., Astron. Astrophys. 447, 31 (2006), astro- 
ph/0510447. 

E. Rasia, S. Ettori, L. Moscaxdini, P. Mazzotta, S. Bor- 
gani, K. Dolag, G. Tormen, L. M. Cheng, and A. Di- 
aferio, Mon. Not. Roy. Astron. Soc. 369, 2013 (2006), 

astro-ph/0602434. 

H. Dahle, N. Kaiser, R. J. Irgens, P. B. Liljc, and S. J. 
Maddox, Astrophys. J. Suppl. 139, 313 (2002). 
R. K. Sheth and G. Tormen, Mon. Not. Roy. Astron. Soc. 
308, 119 (1999), astro-ph/9901122. 

M. S. Warren, K. Abazajian, D. E. Holz, and, L. Teodoro, 

Astrophys. J. 646, 881 (2006), astro-ph/0506395. 

W. L. Freedman et al., Astrophys. J. 553, 47 (2001), 

astro-ph/0012376. 

S. Buries, K. M. NoUett, and M. S. Turner, Phys. Rev. 

D63, 063512 (2001), astro-ph/0008495. 

R. H. Cyburt, Phys. Rev. D70, 023505 (2004), astro- 

ph/0401091. 

P. D. Serpico et al., JCAP 0412, 010 (2004), astro- 

ph/0408076. 

G. Steigman, Int. J. Mod. Phys E15, 1 (2006), astro- 
ph/0511534. 

A. Blanchard, M. Douspis, M. Rowan-Robinson, and 
S. Sarkar, Astron. Astrophys. 412, 35 (2003), astro- 
ph/0304237. 

A. Sandage et al. (2006), astro-ph/0603467. 

D. J. Eisenstein et al., Astrophys. J. 633, 560 (2005), 

astro-ph/0501171. 

A. Lewis and S. Bridle, Phys. Rev. D66, 103511 (2002), 
astro-ph/0205436. 

http://lambda.gsfc.nasa.gov; version v2p2pl. 



